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Abstract: An approximate solution of the Z)-dimensional Schrodinger equa- 
tion with the modified Poschl- Teller potential is obtained with an approxi- 
mation of the centrifugal term. Solution to the corresponding hyper-radial 
equation is given using the conventional Nikiforov-Uvarov method. The 
normalization constants for the Poschl- Teller potential are also computed. 
The expectation values (r~ 2 },(V(r)}, are also obtained using the Feynman- 
Hellmann theorem. 
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1.0 Introduction 

The exact solutions of the Schrodinger equation is only obtainable for 
a few number of physical potentials. However, the Poschl- Teller potential 
happens to be one of the few potentials that enjoy an exact solution of the 
Schrodinger equation. Over the past years, a number of researchers have 
worked on the Poschl- Teller potential [1-5, 11-16, 24]. For instance, a unified 
treatment of the quasi-exactly solvable Poschl- Teller potential was discussed 
by Koc and Koca [15]. The Casimir operator of 50(1,2) and the Poschl- 
Teller potential was presented by Capelas de Oliveira and da Rocha [13] 
using the anti-de sitter(Ads) space time approach. The complex and PT- 
symmetric non-Hermitian Poschl- Teller potential has also been discussed by 
Yesiltas and coworkers [9]. Moreover, the modified Poschl- Teller potential as 
been used to derive the well-known SO (2) spectrum generating algebra for 
an infinite square well problem [44]. Recently, the supersymmetric solution 
of the modified Poschl- Teller was discussed by Aktas and Sever [2], and the 
supersymmetric modified Poschl- Teller and delta-well potential has been 
studied by Diaz et al [5]. 

The extension of physical problems to higher dimensional spaces plays 
an important role in many area of physics. Some of the physical system 
of interest in quantum mechanics which have been thoroughly studied in 
TV-dimensional spaces are the two exactly solvable potentials, which are, the 
Coulomb and harmonic oscillator potentials [28-41]. In particular, Oyewumi 
and Ogunsola [19] have consider the exact solutions of the harmonic oscil- 
lator in multidimensional spaces. Recently, Oyewumi [18] presented the an- 
alytic solutions of the Kratzer-Fues potential in arbitrary dimensions, and 
also the exactly complete solutions of the pseudoharmoinic potential in N- 
dimensions has been discussed by Agboola et al [17]. 

It is therefore the aim of this work to present a multi-dimensional treat- 
ment of the modified Poschl- Teller within the framework of an approxima- 
tion to the centrifugal term. The paper is arranged as follows: the next 
section gives the Schrodinger equation in D-dimensions; section 3 gives a 
brief description of the Nikiforov-Uvarov method, while the following section 
presents the bound-state solutions to the Poschl- Teller potential. In section 
5, some expectation values are calculated; and finallythe conclusions of the 
work are presented in section 6. 
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2. The Schrodinger Equation with the Modified Poschl- Teller 
Potential in /^-dimensions. 

We consider the Schrodinger equation [17-21] 



(1) 



where E is the energy eigenvalue, ^ n ^ m (r, 0) is the wave function and H is 
the Hamiltonian given by 



H = - h - 
2/z 



a-D 9 fn-i 9 



dr 



dr- 



+ V(r), 



(2) 



fj, is the reduced mass, Ti the Planck's constant and V(r) is the modified 
Poschl- Teller potential given as [2,5,11] 



V(r) 



cosh (or) 



(3) 



a being a parameter and Vo is the potential depth. 

Inserting (2)and (3) into (1) and separating the variables as follows 

Eq.(l) reduces to two separate equations 



(4) 



I- 



V 



(k —l)(k — 3) 



cosh 2 (or) J 4:r 2 
with k = D + 2£, and 

A 2 D Y e m (fl) +£(£ + D- 2)Y e m (n) = 
where Y^ m (fi) are the hyperspherical harmonics with 

d 



RnAr) = 0, (5) 



D-2 



a d = e n 



sinf 



(sin 6i) 



i=l \j=o 
and the separation constant 

p = e(e + D-2), 



sinf 



0,1,2,... 



(6) 



d 2 cj) 
(7) 

(8) 



Equation (5) is the hyperradial equation for the modified Poschl- Teller 
potential in D-dimensions, n r the hyperradial quantum number and I is 
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the orbital angular momentum quantum number. 
3. The Nikiforov-Uvarov method. 

In this section, we give a brief description of the conventional Nikiforov- 
Uvarov method. A more detailed description of the method can be ob- 
tained the following references [20] .With an appropriate transformation s = 
s(r), the one dimensional Schrodinger equation can be reduced to a general- 
ized equation of hypergeometric type which can be written as follows: 

y/'( s ) + 44^) + 4r~M s ) = ( 9 ) 

cr[s) <j (S) 

Where cr(s)and a(s) are polynomials, at most second-degree, and f (s)is at 
most a first-order polynomial. To find particular solution of Eq. (9) by 
separation of variables, if one deals with 

^s) = <p(s)y nr ( S ), (10) 

Eq.(9)becomes 
where 



°{s)y'n r + r(s)y' nr + Xy nr = 



i \ ( \ v\°> 

a{s) = 7T{s) wr 

Ha) = f(a) + 2vr(s),r'(s) < 0, 



(11) 
(12) 
(13) 



n(s) = ?L-L ± \/(^) 2 -^+^, (14) 

and 

\ = t + ir'(s). (15) 

The polynomial r(s) with the parameter s and prime factors show the dif- 
ferentials at first degree be negative. However, determination of parameter 
t is the essential point in the calculation of vr(s). It is simply defined by 
setting the discriminate of the square root to zero [20]. Therefore, one gets 
a general quadratic equation for t.The values of t can be used for calculation 
of energy eigenvalues using the following equation 

A = t + 7t'(s) = -n r T>(s) - nr(nr 2 ~ 1} cr"{s). (16) 
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Furthermore, the other part y nr ( s ) of the wave function in Eq.(12) is the 
hypergeometric-type function whose polynomial solutions are given by Ro- 
drigues relation: 

= W)^ {,) ^' )] (17 » 

where B Hr is a normalizing constant and the weight function p(s) must 
satisfy the condition [20] 

(ap)' = rp. (18) 



4.0 Eigenvalues for the modified Poschl- Teller potential. 

For small a, one can use the approximation of the centrifugal potential given 

as 



a 



sinh (or) 



Thus Eq.(5) becomes 
< r e(r) 



2 4 U + 



V 



cosh (or) 



(k - l)(k - 3)q 2 
sinh 2 (ar) 



(19) 



R nr i(r) = (20) 



If we take a transformation s = tanh 2 (ar), after some straight forward 
manipulations, Eq.20 becomes 



1 — 3s 

R nA S ) + —, - R 'nA S ) + 



1 



2s(l-s)" 



4s 2 (l - s) 2 



where 



[-<5 S 2 + ( 7 + 5-e 2 ) S - 7 ] J R nr£ (s) = 

(21) 

(22) 



Comparing Eqs. (21) and (9) we can define the following: 

f(s) = l-3s ,a(s) = 2s(l-s) and a(s) = -5s 2 + (7 + 5 - e 2 )s - 7 (23) 
Inserting these into Eq.(14), we have the following function 

vr(s) = — ^ ± - ^(1 + 45 - 8t)s 2 + (8t - 4(7 + 5 - e 2 ) - 2)s + 4 7 + 1 

(24) 

The constant parameter t can be found by the condition that the discrim- 
inant of the expression under the square root has a double root, i.e., its 
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discriminant is zero. Thus the possible value function for each value of t is 
given as 



tt(s) 



1 - s 



\ [(-2e + yTT4x)s- V / TT4^ for t = -±(7 - 6 + e 2 ) + ±e^TT^y 

\ [(2e + yrq^) s -yrri7l for £ = -±( 7 -<5 + e 2 )-±evrn7~ 

(26) 



By Nikiforov-Uvarov method, we made an appropriate choice of the function 
- \ [(2e + V 1 + 4 7) s - V 71 + 4 7] such that W E q-(!6), we can 



7T(S 



l-s 1 



obtain the eigenvalue equation to be 

-l( n -S+ € 2 )~ € ^T+^~{2e+y/T+^)~ = n r [4+2e+y/T+^]+2n r (n r -l) 



Eq.(27) can be written in the powers of e as follows 

e 2 + e [2(2n r + 1) + + (7 - S) + [(1 + 2n r ) + ^ITi^ 

such that we can obtain the solution as 

' f -2(2n r + l)-yTT4 7 ~+ v / TT45 



(27) 
= 0, (28) 

(29) 



Using the relation in Eq.(22) and k = D + 21, we have the energy eigenvalue 
as 

h 2 a 2 



E„ 



8fi 



2 I 2n r + £- — ) i/l - 



8^0 
/i 2 a 2 



(30) 



One can define the principal quantum number n = 2n r +£, such that Eq.(30) 
becomes 



2^,2 



h a 



(2n + D) + \l + 



h 2 a 2 



,n = 0,l,2,... 



(31) 



We now give the following comments: Eq.(31) is the approximate energy 
spectrum of the Poschl- Teller potential in D-dimensions. However, in order 
to check the validity of the approximation (19), if we consider the s— wave 
and let D = 1, we have the equation 



2^,2 



E n — — 



h a 
~8/7 



n = 0,1,2,... 



(32) 
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which is consistent with the results obtained in previous works [2,24]. Also 
from Eq.(31), one can obtain the critical screening parameter, a c (for which 
E n = 0) as 

/ ~ 
ac ~ Vft 2 [(2n + Df-l] • {66) 

Finally, the curves of the energy levels E n against a for some dimensions 
(figure 1) reveals that the values of the energy for some adjacent dimensions 
becomes closer for higher exited states. 

4.1 Eigenfunctions for the modified Poschl- Teller potential. 



In this section, we obtain the wave functions using the Nikiforov-Uvarov 
method. By substituting ir(s) and a(s) into Eq.(12), and solving the first 
order differential equation to have 



s 2v {\ - s) 



6/2 



2v = 2i + D - 1. 



Also by Eq.(18), the weight function p(s) can be obtained as 

1 



P(s) 



Substituting Eq.(35)into the Rodrigues relation (17), we have 



Vn r {s) = B nr s- V {l - [^ + "(1 - a) 

Therefore, we can write the wave function R nr i(s) as 

RnAs)=Cn r S 2v (l-sy/ 2 P^\l-2s) 



n r +e] 



(34) 



(35) 



(36) 



(37) 



where C„ r is the normalization constant, and we have used the definition of 
the Jacobi polynomials [26], given as 



(-If 



d n 



n!2 n (l - s) a (l + s) b ds n 



(1 



l + s 



(38) 



The curves of the unnormalized hyperradial function are plotted against r 
for some t / states in figures 4-7. 

Moreover, to compute the normalization constant C„ r , it is easy to show 
that 



Jo 



-(g-i) 
r 2 RnAf) 



r D-l dr 



oo rl f ]„ 

Jo 2ay/s(l — s) 

(39) 



= 1 
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where we have also used the substitution s = tanh 2 (ar). Putting Eq.(37) 
into Eq.(39) and using the following definition of the Jacobi polynomial [26] 

Pi a ' b) (s) = ^Vo^ 2Fl a + b + n + l-l + a- 1 -^-), (40) 



we arrived at 



Cn r N nr f 1 s Av ~^{l-sf- 1 [ 2 F l {-n r ,e + v + n r + 1; 1 + v;s)] 2 ds = a (41) 
J o 



r(n r +v+l) 
n r \T(v+l) 



and 2F1 is the hypergeometric function. Using 



where N Hr = \ 

the following series representation of the hypergeometric function 

^(a!,...,^*, ...,<*,) = £ ^Hrr S (42) 

we have 

(43) 

Hence, by the definition of the Beta function, Eq.(43)becomes 



dN n y y { ~ nr)k [ £ + v + n ; ± 1)fc ± ? + n ; ± ^ b ^ + fc + j + 1. 

r (w + l)fefc! + V 2' 

(44) 

Using the relations B(x,y) = ^(al+^f and the Pochhammer symbol 
(a) n = r r(a"' > i Eq.(44) can be written as 

C 2 ^ y (-n r ) fc (e + v + ry + l) fc (4z; + \) k (-n r ) 3 (e + v + ?V + l)j(4^ + k + |)j _ a 
" r " r ^ (e + 4 W + ±) fe (t, + l) fe fc! ^ {e + Av + k + \) 3 {v + l) 3 j\ ~ B{Av + \,e) 

(45) 

Lastly, Eq.(45) can be used to compute the normalization constants for 
n r = 0, 1, 2, ... In particular for the ground state, i.e n r = 0, we have 

c » - (46) 
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5. Some Expectation Values of the Poschl- Teller Potential in 
D-dimensions. 



We now calculate some expectation values of the Poschl- Teller potential 
using the Hellmann-Feynmann theorem (HFT) [27, 28]. Suppose the Hamil- 
tonian H for a particular quantum system is a function of some parameters 
q, and let E(q) and be the eigenvalues and eigenfunctions of H(q) 

respectively, then the HFT states that 



The effective Hamiltonian of the hyperradial function is given as 

V 



H-—f_ ft 2 (21 + P- l)(2l + £>-3) 
2/x dr 2 2fi 4r 2 

In order to calculate (r~ 2 ), we set q = I such that 



cosh (or) 



W)> = |W ^ - 2)<r- 2 > 



and 



dE nA {l) H 2 a 2 



dt 2fi 
Thus by HFT, we have 

v 2 



2 ( 2n r + t + — ) hi/1 f 



h 2 a 2 



a 



2 ( 2n r + £ + — ) + , / I + 



(2£ + D - 2) 

Similarly, putting q = Vq in Eq.(47), we obtain (V(r)} as 

D" 



8fiV 
h 2 a 2 



2 2n r 



+ 2 1+1/1 



?i 2 a 2 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



Finally, we note that the expectation value (T) can also be obtain from the 
fact that E nr = (T) + (V). 
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6. Conclusions. 

In this paper, we present the solutions to the .D-dimensional Schrodinger 
equation with the modified Poschl- Teller potential within the framework 
of an approximation to the centrifugal term. The energy eigenvalues were 
obtained in D-dimensions using the Nikiforov-Uvarov method, and it was 
found to agree with previous works [2, 24]. The corresponding eigenfunc- 
tions of the modified Poschl- Teller Potential were obtained in terms of the 
Jacobi polynomials. The energy levels for some neighboring dimensions were 
found to be getting closer as for higher exited states. Moreover, the graph 
of the absolute value of the unnormalized hyperradial function was plotted 
against r for some I ^ states in 3D and 5D (Figure 2-5), and the curves 
become steeper with increase in the parameter a. 

Finally, the normalization constants were obtained in form of the hyper- 
geometric series and some useful expectation values were computed using 
the Feynman-Hellmann theorem. 
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FIGURES. 

In the figure, U n£ (r) = r-^-^R^r)^ = M = % = 1. 
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Figure 2: \U rir i(r)\ against r for n r = 0,£ = 1 in 3-dimensions. 
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Figure 3: \U rir i(r)\ against r for n r = 0,£ = 1 in 5-dimensions. 
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Figure 5: \U rir i(r)\ against r for n r = 0,£ = 2 in 5-dimensions. 
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